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HIGHLIGHTS 


•  A  battery  signal  pattern  extraction  approach  is  provided. 

•  We  illustrate  that  extracted  patterns  illustrate  the  battery  performance  level. 

•  Clustering  algorithms  are  used  to  demonstrate  visual  interpretations  of  patterns  behaviors. 

•  We  obtain  accurate  shape  classification  performance. 

•  Battery  performances  are  estimated  through  simple  measurements. 
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Direct  estimation  of  battery  performance  is  a  major  challenge  as  ageing  process  is  a  complex  phenom¬ 
enon  not  directly  measurable.  In  this  work  a  new  methodology  is  provided  to  estimate  global  battery 
performances  under  real-life  electric  vehicle  use.  Such  performances  are  estimated  through  battery 
signals  patterns  extraction.  These  signals  patterns  are  used  to  identify  physical  degradation  behavior  of 
batteries. 

The  analysis  framework  is  composed  of  patterns  extraction,  clustering  algorithms,  summarizing  data 
representation  in  the  feature  space  of  cluster  distances  and  classification  algorithms.  This  methodology  is 
then  applied  on  datasets,  acquired  from  batteries  used  on  electric  vehicles,  without  controlled  envi¬ 
ronmental  conditions. 

The  classification  algorithm  accuracy  is  studied  on  the  obtained  real  data.  The  results  suggest  that 
battery  signals  patterns  analysis  provides  an  innovative  technique  for  online  estimation  of  the  battery 
performance  level.  A  detection  of  dysfunctions  caused  by  ageing  is  also  made,  only  based  on  battery 
signals  pattern  extracted  during  real  vehicle  accelerations. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Lithium-ion  (Li-ion)  batteries  are  becoming  the  battery  of  choice 
in  Electric  Vehicles  (EV)  utilization.  However,  battery  health  and 
lifetime  remain  a  major  drawback  to  the  use  of  Li-ion  batteries  in 
stringent  life  requirements.  In  EV  context,  accurate  battery  health 
assessment  is  primordial  to  improve  the  users  confidence  in  the 
battery  range.  Indeed,  it  is  one  of  the  biggest  obstacles  to  wide¬ 
spread  acceptance  of  EVs.  Market  experts  evaluated  the  effects  of 
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low  range  resources  of  EVs,  as  a  significant  feature  for  users'  pur¬ 
chase  intentions  [1], 

The  overall  performance  of  batteries  is  not  constant  along  the 
vehicle  life.  Reduction  of  battery  performances  is  caused  by  various 
internal  and  external  mechanisms  and  is  characterized  by  the  ca¬ 
pacity  fade  as  well  as  an  impedance  augmentation  [2],  For  an 
electric  vehicle  (EV)  utilization,  the  degradation  of  battery  perfor¬ 
mances  can  be  characterized  by  a  diminution  of  the  global  vehicle 
autonomy  available  with  a  full  charge  [3], 

Significant  efforts  have  been  achieved  in  order  to  understand 
the  complex  battery  ageing  process  [4],  The  aim  for  EV  applications 
is  to  estimate  the  performance  level  of  the  battery  with  few  mea¬ 
surements,  in  order  to  obtain  online  estimations  of  the  battery 
performance  level.  Furthermore,  ageing  is  a  complex  phenomenon, 
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difficult  to  estimate  with  only  few  experiments  and  with  online 
constraints.  All  these  constraints  create  an  intricate  compromise 
between  the  model  accuracy  and  its  complexity  [3], 

Different  methods  are  used  to  estimate  the  notion  of  battery 
performance  level.  These  studies  and  methodologies  come  from 
many  various  fields  such  as  electrochemical  modeling  [5,6]  or 
performance  modeling  [7],  The  diversity  and  multitude  of  existing 
studies  dealing  with  battery  ageing  provide  a  large  amount  of  in¬ 
formation.  However,  most  battery  ageing  studies  are  based  on 
direct  factors  dependency  and  have  limited  prediction  ability. 
These  investigations  mainly  rely  on  simulated  data,  under 
controlled  conditions,  which  is  not  totally  representative  of  a  real 
EV  use  [3], 

Moreover,  many  data-driven  methodologies  are  focused  on  the 
battery  capacity  estimation  [8—10],  However,  most  of  these  data- 
driven  approaches  perform  well  on  their  training  data  only,  under 
specific  operational  experiments,  inducing  robustness  and  gener¬ 
alization  mistakes.  In  real  life,  external  conditions  cannot  be 
controlled  and  these  learned  models  are  subject  to  misestimations. 
Thus,  an  accurate  way  of  estimating  battery  capacity  in  real-time 
based  on  real  EV  using  data-driven  algorithm  still  requires  in¬ 
vestigations  [3], 

Battery  performance  changes  can  induce  modifications  of  its 
electrochemical  reactions  in  specific  conditions  11],  as  highlighted 
by  incremental  capacity  analysis  and  differential  voltage  analysis 
[12—14].  Thus,  the  battery  degradations  over  time  lead  to  a  modi¬ 
fication  of  its  behavior  through  the  utilizations  [15,16],  These 
modifications  can  be  detected  with  the  different  battery  measures 
coming  from  real  EV  use.  They  can  provide  indications  about  the 
battery  performance  level.  This  behavior  alteration  is  the  focus  of 
this  study  and  will  be  used  to  estimate  the  battery  performance 
level. 

To  explore  battery  performance  evolution  using  only  real-life 
vehicle  test  data  is  difficult  and  challenging.  The  opportunity 
considered  here  is  to  analyze  the  battery  signals  behavior,  from 
data  collected  from  real-life  EV  operations.  In  our  view,  revisiting 
battery  signals  behaviors  using  pattern  extraction,  analysis  and 
classification  tools  may  provide  new  insights  about  the  rela¬ 
tionship  between  these  signals  and  the  global  battery  perfor¬ 
mance  level,  being  able  to  improve  battery  diagnosis  and 
prognostics. 

In  this  work,  we  propose  an  alternative  approach  by  only  using  a 
data-driven  methodology  developed  from  a  set  of  real  EV  uses. 
Such  a  methodology  requires  a  large  amount  of  training  data  in  the 
development  phase.  In  the  EV  context  this  training  data  require¬ 
ment  is  very  restrictive  and  costly.  To  face  this  problem,  we 
investigate  whether  it  is  possible  to  extract  relevant  features  from 
current  and  voltage  signals  collected  during  real  EV  uses,  under 
uncontrolled  conditions.  A  key  issue  explored  by  this  paper  is  how 
battery  capacity  can  be  estimated  during  real  EV  uses,  without 
specific  requirements,  based  only  on  real  use  data  and  extracted 
features. 

Section  2  presents  the  global  theoretical  framework  and  de¬ 
tails  the  methods  used  for  the  extraction  and  classification  of 
specific  signals  patterns.  This  will  lead  to  an  estimation  of  the 
battery  capacity  level,  during  its  real  uses.  We  study  the  most 
appropriate  extracted  patterns  for  performance  classification,  we 
also  address  the  choice  of  distance  metrics  used  for  comparing 
battery  features.  In  Section  3,  accuracy  of  the  methodology  is 
evaluated  for  several  types  of  classifiers,  using  real  data.  We 
show  that  an  accurate  analysis  of  the  patterns  morphological 
variation,  can  reveal  meaningful  changes  in  the  underlying  bat¬ 
tery  performances.  Finally,  Section  4  presents  a  large  discussion 
of  the  proposed  framework  and  its  applications.  Conclusion  is 
given  in  Section  5. 


2.  Proposed  method 

The  signals  analysis  framework  is  used  in  diverse  application 
fields,  showing  for  example  its  ability  in  classification  tasks  [17,18], 
The  main  applications  of  signal  processing  and  machine  learning 
techniques  are  the  clustering  and  classification  methods,  which  aim 
to  separate  observations  into  groups  in  the  unlabeled  cases  [19], 
Such  algorithms  produce  a  major  step  in  data  analysis  and  explo¬ 
ration.  Cluster  analysis  is  a  common  unsupervised  learning  tech¬ 
nique  for  partitioning  a  dataset  into  subsets  of  data  elements  that 
are  similar  according  to  some  distance  metric  20],  Furthermore, 
based  on  clustering,  groups  analysis  results  led  to  producing 
interesting  interpretations. 

The  pattern  analysis  framework  used  in  this  study  is  illustrated 
in  Fig.  1.  The  first  step  consists  in  a  signal  pattern  extraction  from 
real  measured  battery  data,  to  observe  the  battery  behavior  mod¬ 
ifications.  Then,  different  classification  methods  are  used  to  test 
whether  the  different  signal  morphologies  represent  different 
battery  performance  levels. 

2.3.  Datasets 

Three  EV  instrumented  LiFeCU  (LFP)  batteries  datasets  are  used 
in  this  study.  In  each  dataset,  the  battery  signals  are  acquired  at 
10  Hz  frequency,  during  non-controlled  real  EV  uses,  thereby  data 
representative  of  a  large  variety  of  conditions  an  EV  battery  can  be 
processed.  Moreover,  using  real  data  ensures  compatibility  of 
developed  methodologies  for  embedded  uses.  These  batteries  are 
used  on  three  different  EV,  following  a  vehicle  sharing  concept,  as 
several  drivers  alternatively  take  these  EV  for  their  personal  travels. 
Thus,  these  experimentations  are  representative  of  the  different 
ways  an  EV  can  be  used. 

Table  1  holds  the  main  details  of  the  battery  design  and  their 
characteristics.  Battery  1  and  Battery  3  have  a  final  capacity  lower 
than  Battery  2,  inducing  a  higher  degradation  level  and  thus  a  more 
important  modification  of  their  behavior.  Note  that  the  end  of  life 
(EOL)  criteria  of  the  battery  is  usually  defined  with  a  capacity  level 
equal  to  80%. 

2.2.  Pattern  extraction 

The  methodology  proposed  in  this  work  is  based  on  the 
assumption  of  the  battery  behavior  modification  over  time.  The  aim 
is  to  explore  the  battery  response  depending  on  its  performance 
levels.  For  a  similar  battery  request,  it  is  possible  to  detect  battery 
ageing  effects  based  on  signals  shape  such  as  current  and  voltage 
[12,21],  Thus,  both  current  and  voltage  profiles  present  significant 
behavior  modifications,  throughout  battery  life  requiring  investi¬ 
gation,  characterizing  their  performance  levels. 

However,  these  signal  modifications  are  low  and  complex  to 
identify.  To  observe  a  battery  behavior  alteration,  it  is  necessary  to 
compare  battery  signals  under  comparable  uses.  For  example, 
during  an  identical  speed  profile  criterion,  the  battery  voltage  does 
not  react  the  same  way  depending  on  its  ageing  level. 

To  observe  and  compare  those  signals  modifications,  we  have  to 
consider  only  a  short  length  signal,  with  a  unique  and  common 
dynamic  that  could  be  used  as  a  common  reference.  Hence,  the 
signal  investigation  requires  the  definition  of  criteria  permitting  the 
extraction  of  battery  patterns.  Thus,  duration  and  speed  acceleration 
criteria  define  pattern  extraction  conditions.  These  extraction  criteria 
depend  on  vehicle  and  battery  characteristics  and  have  to  be 
consequently  adapted.  The  impact  of  these  duration  criteria  is  dis¬ 
cussed  further  in  the  following  sections.  Based  on  these  criteria,  the 
obtained  patterns  can  allow  the  identification  of  different  behaviors 
according  to  their  corresponding  battery  performance  levels. 
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Fig.  1.  Battery  signals  analysis  framework  consisting  of  the  data  acquisition,  the  extraction  of  the  patterns,  the  formation  of  a  new  feature  space,  the  patterns  classification  and  the 
estimation  of  the  battery  performance  level. 


The  detection  of  pattern  behavior  alteration  requires  a  long 
battery  solicitation  highlighting  its  performance  level.  Thus,  in  an 
EV  application  the  ideal  conditions  occur  during  a  maximum  ac¬ 
celeration,  allowing  the  signal  patterns  comparison  under  a  com¬ 
mon  reference.  Fig.  2  presents  an  example  of  pattern  extraction 
following  the  criteria  of  an  acceleration  from  20  to  40  km  fr1  in  less 
than  4  s.  Thus  the  pattern  extraction  is  dependent  of  two  criteria: 
speed  thresholds  and  duration  bounds.  Current  and  voltage  pat¬ 
terns  corresponding  to  these  specific  accelerations  are  also 
extracted. 

In  real-life  conditions  it  is  impossible  to  obtain  exactly  twice  the 
same  speed  profile,  implying  a  diversity  among  the  extracted  sig¬ 
nals.  Thus,  the  signal  patterns  are  all  different  and  cannot  be 
directly  compared.  In  order  to  illustrate  the  pattern  variations 
depending  on  the  battery  performance  levels,  average  signals  shape 
are  built  using  the  Filtered  Shape  Averaging  (FSA)  method.  Aver¬ 
aging  results  extracted  from  the  10—60  km  fr1  accelerations  of  the 
Battery  1,  under  four  different  performance  classes  are  presented  in 
Fig.  3.  These  classes  represent  four  disjoint  battery  performance 
levels,  sorted  from  the  less  aged  period  “Class  1"  to  the  more  aged 
battery  level  “Class  4".  Note  that  it  does  not  make  sense  to  compare 
signals  coming  from  different  battery  designs  as  their  behaviors  are 
inherently  very  different. 

The  same  results  for  the  battery  2  and  for  a  shorter  acceleration 
criteria  (20—40  km  fr1)  are  presented  in  Fig.  4. 

Thus,  Figs.  3  and  4  illustrate  variations  of  the  signal  morphol¬ 
ogies  for  different  battery  performance  levels.  For  both  batteries, 
the  extracted  speed  profiles  are  really  close  to  each  other  forming 
good  comparative  samples.  The  corresponding  current  and  voltage 
patterns  are  highly  diverse,  thus  the  four  average  patterns 
demonstrate  a  modification  of  the  battery  behavior,  depending  on 
its  performance  levels.  Indeed,  it  is  clearly  visible  in  Fig.  3  and  in 
Fig.  4  that  the  current  and  voltage  behaviors  are  different  for  each  of 
the  four  performance  classes.  For  example,  the  current  pattern  of 
class  4  (the  most  aged  battery  level)  is  lower  than  the  other  ones, 
but  on  the  contrary  its  corresponding  voltage  pattern  is  at  a  high 
level  at  the  beginning  and  then  decrease  quicker  than  for  the  other 
performance  classes.  These  patterns  also  present  shape  dissimi¬ 
larities  as  visible  in  Fig.  3(b)  illustrating  a  different  shape  behavior 
for  each  of  the  four  classes.  For  example,  average  current  pattern  of 
“Class  1”  decreases  after  6  s,  as  opposed  to  the  growing  average 
current  pattern  associated  to  “Class  4". 

These  behavior  differences  can  be  explained  by  a  slower  reac¬ 
tion  of  the  battery  for  a  same  power  demand  over  its  performances 
level,  but  no  direct  trend  can  be  detected.  However,  interpretations 


Table  1 

Principal  characteristics  of  the  three  studied  batteries. 


Energy 
density 
(Wh  kg-1) 

Nominal 
voltage  (V) 

Total 

capacity 

(Ah) 

Months 
of  use 

Final 

capacity 

(%) 

Battery  1 

80 

400 

75 

59 

87 

Battery  2 

75 

128 

95 

45 

91 

Battery  3 

72 

134.4 

119 

46 

87 

can  be  here  considered  but  are  subject  to  misunderstandings  as  the 
modification  battery  behavior  is  a  complex  problem. 

The  objective  of  the  following  study  is  to  use  these  modifica¬ 
tions  of  the  battery  behavior  to  estimate  the  performance  level, 
only  based  on  the  extracted  patterns. 

2.3.  Dynamic  Time  Warping 

In  order  to  compare  the  signals  pattern  and  quantify  their 
similarities,  we  have  to  consider  a  metric  adapted  to  this  problem. 
Thus,  beyond  usual  measures,  the  current  state-of-the-art  of  shape 
similarity  quantification  is  the  Dynamic  Time  Warping  (DTW).  It 
permits  to  compare  asynchronous  signals  of  different  lengths.  The 
primary  goal  of  DTW  is  to  compare  sequences  respecting  their 
shapes  by  finding  an  optimal  alignment  function  stretching  them. 
Since  its  introduction  in  the  70s,  DTW  has  commonly  been  used  in 
signals  similarity  problems  in  many  fields:  speech  processing,  sig¬ 
nals  recognition,  data  mining  and  imaging  [  22—24], 

This  method  is  based  on  the  Levenshtein  distance  [25]  and  finds 
the  optimal  path  between  two  sequences,  considering  temporal 
distortion.  This  optimal  path  produces  an  alignment  function,  along 
with  a  shape-based  similarity  measure.  Formally,  we  have  two 

sequences  X:  =  (. X\ Xjy )  of  length  N  e  N  and  T:  =  (yi,...,yiw)  of 

length  M  e  N.  In  the  following  we  fix  a  feature  space  denoted  by  T. 
To  compare  two  different  features  x,y  e  F,  one  needs  a  local  cost 
measure,  defined  by  a  function  c: 

C  :  J~  X  J~  —*  R>o  ( 1 ) 

Typically,  the  cost  c(x,y)  is  low  if  x  and  y  are  similar  to  each  other, 
otherwise  c(x,y)  is  high.  Evaluating  the  local  cost  measure  c(x,y)  for 
each  pair  of  elements  of  the  sequences  X  and  V,  one  obtains  the  cost 
matrix  C  e  rWxM  defined  by  C(ij)  =  c(Xj,yj).  The  goal  is  to  find 
the  alignment  between  X  and  Y  minimizing  the  overall  cost. 
A  warping  path  is  a  sequence  p  =  (P\,.-,Pl)  with 

Pi  =  (n(,  m()e[  1  :  N]  x  [1  :  M],  V/e[l  :  L\,  satisfying  the  following 
conditions: 


0  10  20  30  40  50 


Time  (minutes) 

Fig.  2.  Patterns  extraction  of  an  acceleration  from  20  to  40  km  h  \  from  EV  battery  2 
speed  signal,  corresponding  current  and  voltage  patterns  are  also  extracted. 
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Fig.  3.  (a).  Speed  profile  -  average  patterns  profile  corresponding  to  a  maximum  acceleration  under  four  different  performance  levels  of  battery  1,  based  on  10-60  km  h_1  ac¬ 
celerations.  (b).  Current  profile  -  average  patterns  profile  corresponding  to  a  maximum  acceleration  under  four  different  performance  levels  of  battery  1,  based  on  10-60  km  h_1 
accelerations,  (c).  Voltage  profile  -  average  patterns  profile  corresponding  to  a  maximum  acceleration  under  four  different  performance  levels  of  battery  1,  based  on  10-60  km  h-1 
accelerations. 


r  Pt  =  (1,1)  and  pL  =  (N,M) 

<  nj  <  ...  <  rii  and  irq  <  ...  <  mL  (2) 

1  Pf+i  -  Pie{(l>0),  (0, 1),  (1, 1)},  Vle[l  :  L  —  1] 

A  warping  path  p  =  (pi,.  ,.,pi)  defines  an  alignment  between  two 
sequences  X  and  Y  by  assigning  the  element  xn,  of  X  to  the  element 
ym,  of  Y.  The  alignment  conditions  imply  that  the  first  elements  of  X 
and  Y  as  well  as  their  last  elements  are  aligned  to  each  other.  The 
total  cost  Cp(X,Y)  of  a  warping  path  p  between  X  and  Y  with  respect 
to  the  local  cost  measure  c  is  defined  as: 


L 

Cp(X,Y)  =  ^c(xn„ym,)  (3) 

;=i 

Examples  of  paths  between  two  signals  X  =  {xp..  ,,xg}  and  Y= 
{yi,...,y8}  are  given  in  Fig.  5,  illustrating  an  admissible  path 
following  the  Equation  (2)  (Fig.  5a)  and  two  non  admissible  paths 
due  to  boundary  conditions  (Fig.  5b)  and  to  step  size  conditions 

(Fig.  5c). 

Furthermore,  an  optimal  warping  path  between  X  and  Y  is  a 
warping  path  p*  minimizing  total  cost  among  all  possible  warping 


Fig.  4.  (a).  Speed  profile  -  average  patterns  profile  corresponding  to  a  maximum  acceleration  under  four  different  performance  levels  of  battery  2,  based  on  20-40  km  h  1  ac¬ 
celerations.  (b).  Current  profile  -  average  patterns  profile  corresponding  to  a  maximum  acceleration  under  four  different  performance  levels  of  battery  2,  based  on  20—40  km  h  1 
accelerations,  (c).  Voltage  profile  -  average  patterns  profile  corresponding  to  a  maximum  acceleration  under  four  different  performance  levels  of  battery  2,  based  on  20-40  km  h  1 
accelerations. 
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*i  x2  x3  x4  x5  x6 

Signal  X 


*1  *2  X3  X4  X5  X6 
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Fig.  5.  (a).  Admissible  warping  path  —  illustration  of  paths  of  index  pairs  for  a 
sequence  X  =  {xi,...,X6}  and  a  sequence  Y  =  {yi,...^}-  (b).  Example  of  a  non  admissible 
warping  path  due  to  boundary  condition  -  illustration  of  paths  of  index  pairs  for  a 
sequence  X  =  {xi,...,x6}  and  a  sequence  Y  =  (c)-  Example  of  a  non  admissible 

warping  path  due  to  step  size  condition  —  illustration  of  paths  of  index  pairs  for  a 
sequence  X  =  {xi,...,X6}  and  a  sequence  Y  =  {yi.-.-O's}* 


paths.  The  DTW  distance  doTw(X^)  between  X  and  Y  is  then  defined 
as  the  total  cost  of  the  optimal  warping  path  p*: 

dDTw(X,  Y)  =  cp.  (X,  Y)  =  min{cp(X,  Y) |  Vp}  (4) 

This  DTW  distance  permits  the  comparison  and  the  quantifi¬ 
cation  of  different  signals  shape,  it  is  particularly  adapted  to  battery 
signals  evolution.  Therefore,  this  distance  measures  the  difference 
between  each  extracted  pattern. 


2.4.  Classification  methods 


Based  on  the  DTW  distance  applied  on  the  extracted  patterns, 
the  aim  is  to  estimate  the  battery  performance  level.  In  this  section, 
different  classification  and  clustering  methods  are  presented.  All  of 
these  methods  require  a  distance  permitting  the  comparison  be¬ 
tween  observations,  which  is  going  to  be  the  DTW  distance  in  our 
application  as  it  is  adapted  to  the  nature  of  the  analyzed  patterns. 

The  input  of  the  clustering-classification  framework  is  a  dataset 
B  =  {(f>i,/i),...,(t>n,fn)}.  composed  of  n  signals  pattern  bi  along  with 
their  associated  performance  label  1,-. 

The  considered  dataset  B  is  partitioned  into  a  training  set  Btraln 
and  a  testing  set  Btest,  permitting  to  build  classification  algorithms 
based  on  the  training  dataset.  The  following  detailed  classification 
algorithm  provides  a  label  estimation  /<■  for  an  observation  b;  e  Btest. 
Since  the  class  labels  represent  the  battery  performance  level,  the 
classification  accuracy  CA  on  a  method  is  defined  as  the  percentage 
of  data  elements  classified  in  their  actual  label  /,■: 


r  a  _  LL _ ! _ L 

|Btest| 


(5) 


2.4. 1.  K-Nearest  Neighbors 

The  X-Nearest  Neighbors  (X-NN)  approach  has  been  widely 
investigated  and  has  shown  to  be  a  powerful  non-parametric 
technique  in  classification  problems.  This  algorithm  is 


computationally  simple  and  commonly  used  as  a  reference,  thus  a 
new  object  is  assigned  to  the  majority  class  voted  among  its  X 
nearest  neighbors  in  the  training  set  [26].  Formally,  given  a  labeled 
training  set  Btraln  and  a  test  pattern  b;  c  Btest,  X-NN  classifies  by 
choosing  X  train  data  elements  {bi,. . .,b/<-}  c  Btraln  that  are  the  X 
closest  neighbors  of  b;  under  the  used  metric.  Then,  b ,•  will  be 
classified  as  the  statistical  mode  from  the  labels  of  the  X  nearest 
neighbors  {Zlt. . .,//<-}-  Note  that  the  parameter  X  is  a  positive  integer 
strongly  influencing  the  algorithm  accuracy. 

2.4.2.  K-Means 

The  X-Means  algorithm  is  one  of  the  most  popular  and  widely 
used  as  a  clustering  method.  This  algorithm  follows  a  local  opti¬ 
mization  strategy.  It  starts  with  an  initial  partition  of  the  n  obser¬ 
vations  of  Btiain  into  X  clusters,  and  then  uses  an  iterative  process  to 
optimize  the  clusters  quality. 

The  aim  is  to  group  n  observations  into  X  separate  clusters 
{Ci,...,C«j,  so  that  the  within-cluster  sum  of  squares  is  minimized. 
This  distortion  E  is  measured  by  E  =  ]T)"=1  (  min  D(b;,  Q)),  with  D 
the  DTW  distance  measure.  The  distance  feelwfcn  an  element  b, 
and  a  cluster  Ck  is  defined  as  the  distance  between  bi  and  the  center 
of  this  cluster  cf  at  the  step  t.  This  center  of  a  cluster  is  built  ac¬ 
cording  to  the  Filtered  Shape  Averaging  (FSA)  method,  producing 
an  average  signal  representative  of  the  class  behavior,  using  the 
DTW  optimal  path  with  an  iterative  process.  Thus,  this  methodol¬ 
ogy  is  dependent  on  the  choice  of  initial  positions  of  cluster  centers 
(c^°\  ...,c^).  In  this  study  these  initial  positions  are  randomly 
chosen  among  the  observations.  At  the  end  of  the  iterative  proce¬ 
dure,  all  observations  of  the  training  set  Btraln  are  assigned  to  a 
single  cluster. 

With  the  formation  of  the  X  clusters,  a  new  observation  b;  e  Btest 
is  associated  to  the  nearest  cluster  Cj.  Thus,  the  label  associated  to 
the  pattern  b;  is  the  statistical  mode  from  the  labels  of  this  cluster 

Cj- 

2.4.3.  Hierarchical 

Hierarchical  clustering  is  applied  to  the  n  training  patterns  of 
Btram  t0  Qbtgjjj  a  partition  of  this  dataset  into  disjoint  clusters,  in  a 
way  that  data  elements  within  the  same  cluster  share  a  similarity 
according  to  the  chosen  DTW  distance.  This  clustering  is  done 
following  an  agglomerative  hierarchical  clustering  procedure  [27] 
that  initially  partitions  a  set  of  n  data  elements  into  X  distinct 
clusters,  each  containing  one  data  element,  and  then  iteratively 
merges  the  two  most  similar  clusters,  until  the  entire  dataset  is 
aggregated  into  a  unique  cluster.  The  created  hierarchical  tree  can 
next  be  pruned  so  that  the  number  of  clusters  X  is  obtained. 

Thus,  data  elements  below  each  cut  are  assigned  to  a  single 
cluster,  creating  the  output  data  partitioning  to  clusters  {ci,...,ck}. 
The  algorithm  requires  a  cluster  similarity  criterion  for  choosing 
the  next  two  clusters  to  be  merged.  We  have  simply  used  the 
average  similarity  between  all  the  centers  of  both  clusters. 

Following  the  same  classification  process  as  the  X-Means  algo¬ 
rithm,  a  new  observation  is  associated  to  the  label  of  its  nearest 
cluster. 

2.4.4.  Support  Vector  Machine 

The  theoretical  foundations  of  Support  Vector  Machine  (SVM) 
were  developed  by  Vapnik  [28]  the  method  has  been  highly 
explored  during  recent  years.  The  SVM  method  is  applied  for 
classification  in  various  domains  providing  excellent  results 
[29—32].  Using  the  set  of  training  patterns  Btraln,  and  their  corre¬ 
sponding  class  identifiers;  it  assigns  the  correct  class  to  a  newly 
presented  object.  For  the  task  of  learning,  a  measure  of  similarity 
between  the  objects  of  Btialn  is  necessary,  so  that  classifiers  are 
created  separating  the  data  space  into  separate  areas.  A  reasonable 
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measure  of  similarity  is  of  the  form  k  :  Btraln  x  Btraln— >R,  where  k  is 
called  a  kernel.  In  our  application,  the  chosen  kernel  k  is  the  DTW 
Kernel  (DTWK)  [33],  permitting  the  use  of  the  DTW  metric  in  the 
SVM  method. 

Note  that  the  SVM  classification  is  initially  adapted  to  a  two- 
class  problem.  In  order  to  adapt  the  SVM  to  a  multi-class  we 
reduce  the  single  multiclass  problem  into  multiple  binary  classifi¬ 
cation  problems  [34], 

3.  Results  on  real  data 

The  described  framework  is  applied  separately  on  each  battery 
dataset.  Several  extraction  criteria  are  studied,  permitting  the 
comparison  of  the  results  depending  on  the  patterns.  Thus,  the 
global  methodology  is  applied  on  two  types  of  speed  profile,  one 
considering  an  acceleration  from  20  to  40  km  hr1,  and  the  other 
from  10  to  60  km  hr1.  Notice  that  the  longest  acceleration  contains 
more  information.  However,  longer  acceleration  profiles  require 
larger  datasets  to  obtain  enough  patterns  for  the  methodology 
process.  These  extracted  patterns  are  then  used  in  the  global  al¬ 
gorithm  and  the  results  are  explored. 

3.1.  Clustering 

The  first  step  consists  in  the  exploration  of  the  patterns  clus¬ 
tering  results.  Hierarchical  clustering  and  K-Means  methodologies 
can  be  applied  to  the  extracted  patterns  described  above.  Cluster 
analysis  is  able  to  identify  distinct  patterns  morphologies.  Fig.  6 
presents  for  example  a  clustering  based  on  current  patterns 
extracted  during  specific  accelerations,  with  the  average  current 
patterns  of  each  cluster  resulting  of  the  FSA  methodology.  These 
results  illustrate  the  diversity  of  the  battery  current  shape  between 
the  various  8  clusters.  For  example,  cluster  8  presents  decreasing 
currents  after  7  s,  as  opposed  to  cluster  4  containing  constantly 
growing  current  patterns.  Note  that  the  speed  profile  correspond¬ 
ing  to  the  current  patterns,  as  previously  explained,  does  not  pre¬ 
sent  significant  differences  due  to  the  pattern  extraction  mode.  This 
result  permits  to  demonstrate  the  variability  of  the  current  patterns 
behavior. 

As  the  data  are  provided  from  real-life  tests,  these  differences 
can  be  caused  by  external  conditions  such  as  temperature,  State  Of 
Charge  (SOC)  battery  level.  Indeed,  a  cluster  can  for  example  be 


composed  of  pattern  measured  during  summer  only,  whatever  the 
battery  performance  level. 

To  demonstrate  the  independence  between  the  cluster  compo¬ 
sition  and  the  external  conditions,  Fig.  7(a)  and  (b)  presents  the 
observations  distribution  over  temperature  and  SOC  conditions, 
with  their  associated  cluster.  Fig.  7  illustrates  the  homogeneous 
distribution  of  temperature  and  SOC  over  clusters.  Thus,  given  a 
temperature  or  a  SOC  it  is  impossible  to  predict  its  associated 
cluster.  To  quantify  these  relationships,  the  absolute  correlation 
coefficient  between  the  temperature  and  the  selected  cluster 
mainly  lower  than  0.1,  as  well  as  the  correlation  between  SOC  and 
clusters.  These  results  clearly  illustrate  the  global  variability  of 
external  conditions  in  all  clusters,  and  the  independence  between  a 
cluster  composition  and  their  corresponding  external  conditions. 
We  can  conclude  that  external  conditions  have  no  influence  on  the 
results. 

Based  on  these  cluster  investigations  the  battery  behavior  on  a 
small  acceleration  part  permits  the  identification  of  patterns  shape 
differences  according  to  its  performance  level  regardless  of  the 
external  conditions. 

It  is  also  considered  to  use  this  approach  to  identify  behavior 
characteristics  of  a  battery  fault.  This  fault  detection  can  be  made  by 
a  simple  outlier  identification  with  the  use  of  this  clustering  work. 
Thus,  a  battery  default  will  be  clearly  visible  as  the  corresponding 
signal  patterns  will  be  isolated  in  a  unique  cluster.  This  method¬ 
ology  consequently  permits  a  rapid  diagnosis  of  the  global  battery 
behavior  to  detect  abrupt  changes  of  battery  behavior  representa¬ 
tive  of  a  malfunction. 

3.2.  Classification 

The  clustering  investigation  induces  a  good  separation  of  the 
patterns  studied,  according  to  their  respective  level  of  performance. 
The  idea  is  now  to  use  these  modifications  of  the  battery  behavior 
to  predict  its  level  of  performance,  only  based  on  a  pattern 
extracted  from  a  single  acceleration.  The  four  classification 
methods  described  in  Section  2.4  (K-NN,  K-Means,  Hierarchical  and 
SVM)  are  used  and  their  results  are  compared  based  on  the 
different  datasets. 

Three  datasets  coming  from  three  different  batteries  are 
explored,  as  explained  in  Section  2.1.  Each  dataset  corresponds  to  a 
real  EV  use.  In  order  to  compare  the  classification  methods 


Cluster  1  (26)  Cluster  2  (13)  Cluster  3  (24)  Cluster  4  (18)  Cluster  5  (12)  Cluster  6  (10)  Cluster  7  (6)  Cluster  8  (5) 
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Fig.  6.  Clustering  result  of  the  K-Means  method  for  10-60  km  h  1  accelerations  of  battery  1  -  the  number  of  patterns  in  each  cluster  are  in  parenthesis,  in  red  is  presented  an 
average  current  pattern  of  each  cluster,  obtained  by  FSA.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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Fig.  7.  (a).  Temperature  repartition  of  the  clusters  -  distribution  of  the  external  conditions  of  each  observation,  along  with  their  selected  cluster  based  on  a  /C- Means  clustering  with 
the  battery  1.  (b).  SOC  repartition  of  the  clusters  —  distribution  of  the  external  conditions  of  each  observation,  along  with  their  selected  cluster  based  on  a  /(-Means  clustering  with 
the  battery  1. 


accuracy,  the  methodology  is  applied  to  different  sets  of  extracted 
patterns  corresponding  to  a  20—40  km  fr1  or  to  a  10—60  km  fr1 
acceleration.  Table  2  presents  the  number  of  extracted  patterns  by 
battery  during  their  entire  experimentations,  for  two  different  ac¬ 
celeration  profiles  criteria.  For  example,  on  Battery  1, 143  patterns 
satisfy  the  extraction  criteria  of  a  10—60  km  hr1  acceleration  during 
a  specific  duration,  depending  on  the  vehicle  characteristics, 
defined  by  experts.  Note  that  the  number  of  extracted  patterns  is 
highly  dependent  of  the  battery  uses,  explaining  the  number  of 
patterns  corresponding  to  20—40  km  h  1  accelerations  in  Battery 
experiments.  Thus,  extraction  process  provides  a  large  amount  of 
patterns  corresponding  to  an  identical  speed  profile. 

Each  dataset  B  =  {(fq,/i),... ,(&„,/„)}  is  composed  with  different 
extracted  patterns  i>;,  associated  to  a  class  of  battery  performance 
These  datasets  B  are  then  separated  into  a  training  part  Btram  con¬ 
sisting  of  75%  of  the  entire  dataset  B,  and  into  a  test  set  Btest 
composed  with  the  other  25%  of  the  dataset  B. 

As  explained  in  Section  2.4,  each  classification  method  builds  a 
model  based  on  a  training  part  Btraln  of  this  dataset.  The  learned 
model  is  then  applied  to  the  test  patterns  b ;  e  Btest,  and  the  clas¬ 
sification  accuracy  CA  is  the  proportion  of  well  classified  test  ele¬ 
ments.  Note  that  for  each  method,  the  parameters  are  defined 
following  a  cross  validation  step  using  the  training  dataset.  Hence, 
the  parameters  are  automatically  defined  during  the  preprocessing 
step. 

In  this  section  we  apply  classification  methodologies  to  extrac¬ 
ted  voltage  and  current  patterns,  and  to  their  corresponding  power 
pattern. 

3.3.  Two  classes 

Consider  here  a  binary  classification  for  extracted  patterns.  For 
each  battery,  the  two  disjoint  classes  correspond  to  a  performance 
class  level.  Thus,  one  class  contains  the  observations  associated 
with  the  first  quarter  of  the  performance  level  and  the  other  class 
holds  a  last  quarter  of  battery  performances.  It  means  that  the  first 
quarter  of  the  battery  life  is  opposed  to  the  last  quarter,  each  of 
them  composing  a  specific  performance  class.  Thus,  for  a  battery, 
this  classification  uses  only  half  of  its  total  extracted  patterns. 
Considering  a  test  pattern,  its  corresponding  battery  performance 


Table  2 

Number  of  extracted  patterns  for  each  battery  with  two  different  acceleration 
criteria. 


10-60  kmh1 

20-40  km  h  1 

Battery  1 

143 

5192 

Battery  2 

139 

662 

Battery  3 

249 
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class  is  estimated  by  each  of  the  presented  classification  methods. 
Table  3  presents  a  comparison  of  the  accuracies  of  each  classifica¬ 
tion  method,  for  the  three  batteries  and  for  two  different  acceler¬ 
ation  criteria. 

The  presented  results  of  Table  3  clearly  demonstrate  the  ability 
of  the  methodology  to  perform  with  a  good  accuracy  in  the  battery 
class  of  performance  identification.  Indeed,  the  best  accuracy  is 
obtained  with  the  10—60  km  fr1.  This  can  be  explained  because  as 
the  signal  gets  larger  it  contains  more  and  more  information  on  the 
battery  behavior.  Thus,  the  best  performance  is  obtained  with 
Battery  3  over  10—60  km  lr1  accelerations  with  the  1C-NN  method, 
based  only  a  training  dataset  composed  of  less  than  200  patterns. 

These  results  also  point  out  the  AT-NN  and  the  SVM  as  the  most 
accurate  methods.  These  methods  yield  a  correct  battery  perfor¬ 
mance  class  for  a  new  signal  pattern  of  around  75%,  depending  on 
the  patterns  lengths  and  the  dataset  size.  Note  that  the  classifica¬ 
tion  accuracy  is  up  to  90%  for  the  /C-NN  method  based  on  voltage 
patterns  extracted  from  10—60  km  h  This  result  is  presented  is 
Fig.  8,  illustrating  the  methodology  accuracy.  Compared  to  a 
random  process  giving  a  50%  accuracy  by  assigning  to  each  pattern 
one  of  the  two  classes  randomly,  the  presented  methodology 
demonstrates  its  ability  to  estimate  battery  performance  level. 

Note  that  the  performances  obtained  from  Battery  2  are  lower 
than  other  batteries,  which  is  explained  by  the  fact  that  this  battery 
has  a  higher  final  capacity  than  the  two  others  batteries,  providing 
fewer  patterns  due  to  the  low  utilizations  of  this  Battery  2  and 
consequently  a  reduced  set  of  training  data  Btraln.  Moreover,  this 
Battery  2  was  less  aged  than  the  others  batteries,  inducing  a  smaller 
modifications  of  its  behavior.  The  signals  alterations  are  more 
complex  to  detect  than  for  a  battery  with  a  long  utilization  life. 
Consequently  its  behavior  modifications  are  less  sensible  and  the 
patterns  alterations  are  more  difficult  to  detect. 


Table  3 

Classification  accuracy  of  a  two-class  classification  of  the  battery  performance  level 
for  each  pattern  (%)  —  in  bold  the  best  accuracies  (>75%). 


Bat.  1 

Bat.  2 

Bat.  3 

I 

U 

P 

I 

U 

P 

I 

U 

P 

20-40  km  h~] 
K-  NN 

72 

77 

74 

65 

60 

56 

69 

89 

69 

K-  Means 

67 

64 

66 

63 

62 

60 

69 

81 

69 

Hierarchical 

61 

63 

60 

56 

62 

60 

67 

85 

69 

SVM 

69 

70 

71 

64 

62 

61 

67 

86 

67 

10-60  km  h 

K-NN 

79 

59 

69 

66 

80 

69 

80 

90 

78 

K-  Means 

79 

72 

76 

60 

71 

71 

71 

75 

71 

Hierarchical 

60 

62 

59 

54 

69 

69 

69 

77 

71 

SVM 

72 

72 

69 

63 

69 

63 

76 

80 

75 
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Fig.  8.  Prediction  distribution  of  the  performance  class  by  the  K-NN  method,  based  on 
voltage  patterns  of  10-60  km  h  1  acceleration  of  battery  3. 

Using  the  K-NN  approach  applied  to  Battery  3,  we  propose  a 
sensitivity  study  of  the  pattern  extraction  criteria.  We  consider  the 
fix  speed  accelerations  thresholds  at  10—60  km-1,  and  we  use 
diverse  duration  bounds  to  extract  corresponding  voltage  patterns. 
Thus  for  each  acceleration  duration  bounds,  the  K-NN  methodology 
creates  a  model  based  on  the  extracted  patterns.  All  of  these  tests 
allow  the  production  of  an  accuracy  level  of  the  formed  model. 
Fig.  9  illustrates  the  model  accuracy  evolution  for  different  accel¬ 
eration  duration  bounds. 

Fig.  9  illustrates  the  compromise  in  the  choice  of  extraction 
criteria.  We  can  see  the  possibility  to  improve  the  accuracy  pre¬ 
sented  in  Table  3  from  90%  to  a  98%  accuracy,  by  considering 
reduced  acceleration  duration  bounds.  This  bound  definition  is  an 
important  part  of  the  pattern  extraction  methodology  as  the 
extraction  criteria  are  currently  defined  a-priori  from  expert 
advices. 

3.4.  Four  classes 

With  a  good  global  accuracy  in  a  two-classes  classification,  the 
aim  is  to  use  the  same  methodology  to  predict  with  a  greater 
precision  the  battery  performance  level.  Thus,  Table  4  presents  the 
classification  rates  of  4-class  classifications.  The  interpretation  of 
these  results  is  similar  to  the  conclusion  made  from  the  2-class 
results.  Indeed,  the  best  classification  performances  are  obtained 
with  the  K-NN  method,  for  the  10—60  km  fr1  accelerations.  Note 
that  a  random  classification  gives  here  an  accuracy  of  25%,  thereby 
affirming  that  the  global  methodology  permits  the  estimation  of 
the  battery  performance  level  with  a  good  success  rate.  It  is 
important  to  observe  the  low  results  obtained  by  the  SVM  methods, 
opposed  to  the  performances  of  the  previous  2-class  approach.  This 
comment  is  here  explained  by  the  definition  of  the  used  kernel  in 
the  SVM  process.  Indeed  the  defined  DTW  Kernel  (DTWK)  is  not 
positive  definite  in  theory,  but  is  commonly  p.d.  in  practical  ex¬ 
periments.  However,  the  results  presented  in  Table  4  demonstrate 
the  limitation  of  this  kernel's  utilization. 


Fig.  9.  Evolution  of  the  two-classes  K-NN  accuracy  from  patterns  extracted  from 
10-60  km  h  1  accelerations,  of  battery  3,  with  different  duration  bounds. 


Table  4 

Classification  accuracy  of  a  four-class  classification  of  the  battery  level  of  perfor¬ 
mance  for  each  datasets  (%)  —  in  bold  the  best  accuracies  (>50%). 


Bat.  1 

Bat.  2 

Bat.  3 

1 

U 

P 

I 

U 

P 

I 

U 

P 

20-40  km  Tr1 
K-NN 

47 

42 

48 

29 

35 

33 

41 

62 

41 

K-  Means 

40 

37 

41 

35 

42 

33 

49 

69 

42 

Hierarchical 

31 

30 

32 

30 

31 

31 

45 

58 

44 

SVM 

37 

36 

37 

33 

32 

30 

43 

55 

40 

70-60  km  h-1 
K-NN 

48 

38 

48 

46 

49 

43 

51 

67 

48 

/(-Means 

48 

48 

52 

49 

46 

40 

44 

55 

42 

Hierarchical 

45 

31 

35 

28 

40 

40 

41 

57 

34 

SVM 

48 

34 

38 

37 

46 

34 

32 

42 

33 

The  results  of  a  4-class  K-NN,  based  on  voltage  patterns  from  the 
10—60  km  fr1  of  Battery  3,  are  presented  in  Fig.  10.  This  Fig.  10 
illustrates  the  ability  of  the  proposed  methodology  to  estimate 
most  of  times  the  correct  battery  performance  class,  only  based 
here  on  the  extracted  voltage  pattern.  For  example,  the  observa¬ 
tions  corresponding  to  the  battery  performance  class  1  are  strongly 
well  classified.  However,  we  can  observe  mispredictions,  especially 
between  close  classes,  but  rarely  for  remote  classes. 

This  representation  demonstrates  the  independence  between 
the  prediction  made  and  the  unused  variables.  This  independence 
of  clustering  on  battey  temperature  and  SOC  is  illustrated  by  the 
homogeneous  distribution  of  each  of  these  conditions  in  all  four 
composed  classes.  For  example,  the  temperature  distribution  con¬ 
tained  into  the  class  2  contains  sample  associated  with  low  values 
(-0  °C)  as  well  as  high  temperatures  (-40  °C).  That  is  to  say  that  the 
built  model  equally  performs  for  all  values  of  the  State  Of  Charge 
(SOC)  and  of  the  internal  battery  temperature. 

Thus,  the  presented  results  demonstrate  the  good  precision  of 
the  methodology  as  a  unique  voltage  pattern  permits  to  estimate 
the  battery  performance  level  with  an  interesting  accuracy. 
Furthermore  it  is  interesting  to  observe  the  different  accuracies 
depending  on  the  used  signal.  Indeed,  the  voltage  patterns  provide 
the  best  results,  inducing  that  the  battery  behavior  modifications 
are  more  detectable  inside  the  voltage  signals  than  in  current  sig¬ 
nals.  On  the  contrary,  the  power  patterns  lead  to  lower  accuracy 
rate.  This  point  can  be  explained  as  the  power  is  defined  as  the 
product  between  current  and  voltage.  Thus  both  behaviors  of  cur¬ 
rent  and  voltage  can  be  attenuated  or  altered  by  this  product.  For 


Fig.  10.  Prediction  repartition  of  the  performance  class  by  the  K-NN  method  with 
battery  3,  depending  on  the  temperature  and  SOC  values. 
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the  best  estimation  precision,  the  presented  results  tend  to  advo¬ 
cate  the  choice  of  voltage  patterns  in  the  global  methodology. 

In  order  to  explore  the  methodology  performances  on  smaller 
performances  classes,  Fig.  11  presents  the  results  of  the  K-NN 
methods  on  an  8-class  classification,  only  based  on  voltage  patterns 
from  20—40  km  h  .  The  choice  of  these  acceleration  criteria  is 
made  to  obtain  enough  training  patterns  to  perform  in  a  8-classes 
approach. 

This  result  demonstrates  the  good  detection  of  battery  behavior 
modifications  depending  on  its  performance  level,  only  based  on 
patterns  of  voltage  extracted  on  a  simple  acceleration.  Indeed,  most 
class  estimations  are  correct  or  are  a  neighbor  of  the  real  class.  For 
example,  in  Fig.  11  misclassed  Class  8  patterns  are  associated  to 
Class  7,  corresponding  to  the  nearest  performance  class  of  Class  8. 
The  difference  between  the  battery  performance  classes  is  here  less 
than  2%  of  the  battery  capacity,  inducing  the  good  identification  of 
the  battery  behavior  alteration. 

4.  Discussion 

The  proposed  statistical  approach  provides  an  innovative 
framework  to  use  patterns  from  battery  signals.  This  study  does  not 
consider  physical  hypothesis  allowing  an  adaptation  of  the  global 
process  with  any  battery  technology  and  perform  in  real  uses 
without  specifications,  demonstrating  the  powerful  of  signal  pat¬ 
terns  investigation. 

Simple  patterns  extraction  and  their  analysis,  presented  in 
Section  2.2  permit  the  identification  of  a  battery  behavior  modifi¬ 
cation  over  different  performance  levels  (Figs.  3  and  4).  Then, 
deeper  signal  patterns  investigation,  introduced  in  Section  2.4, 
demonstrates  the  ability  of  the  patterns  extraction  to  highlight  the 
battery  behavior,  only  based  on  the  signals  coming  from  real  ac¬ 
celerations.  Indeed,  the  results  of  Section  3  confirm  the  interest  of 
this  methodology  in  the  estimation  of  the  battery  performance 
level,  and  for  fault  detection. 

The  four-classes  classification  results  presented  in  Table  4 
illustrate  the  good  accuracy  of  the  methodology  to  detect  battery 
behavior  over  the  performance  level,  only  based  on  the  battery 
signal  patterns  measured  during  vehicle  uses,  without  requirement 
or  specification.  The  obtained  results  present  diverse  accuracy  for 
the  selected  patterns  extraction  mode,  thus  it  is  of  paramount 
importance  to  define  optimal  patterns  extraction  criteria,  among  all 
accelerations  profiles,  illustrated  by  Fig.  9.  This  study  highlights  the 
importance  of  the  pattern  extraction  criteria  definition,  as  they 
directly  impact  the  accuracy  of  the  developed  methodology.  In  this 
presented  experiment  all  these  extractions  criteria  are  preliminary 
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Fig.  11.  Estimation  repartition  of  the  8-classes  K-NN  method,  based  on  voltage  patterns 
of  20-40  km  h  1  acceleration  of  battery  3  -  with  a  55%  accuracy. 


defined  manually,  however  Fig.  9  demonstrates  the  significance  of 
exploring  other  definition  modes  to  detect  optimal  criteria. 

Table  2  details  the  number  of  extracted  patterns  for  two 
different  acceleration  profiles,  illustrating  the  diversity  of  extracted 
patterns  according  to  speeds  considered.  Indeed,  a  10—60  km  It1 
acceleration  pattern  is  harder  to  obtain  during  real  EV  uses  than  a 
20—40  km  tr1  speed  profile. 

Thereby,  the  presented  experiments  provide  result  illustrations 
with  patterns  extracted  from  the  10—60  km  h  1  accelerations. 
Indeed,  longer  patterns  contain  more  battery  behavior  information 
and  thus  permit  a  better  identification  of  the  battery  performance. 
This  compromise  between  frequency  and  accuracy  of  the  estima¬ 
tions  has  to  be  selected  in  accordance  with  the  applications  con¬ 
straints.  Consequently,  better  estimation  results  are  obtained  with 
Battery  3  containing  the  most  extracted  patterns,  with 
10—60  km  hr1  accelerations  criteria.  Among  the  detailed  classifi¬ 
cation  methods,  the  K-NN  approach  is  associated  to  the  best  esti¬ 
mations  accuracies  inducing  a  preference  for  the  application  of  this 
method.  However,  the  SVM  algorithm  permits  to  obtain  high  ac¬ 
curacies  in  a  two-class  context  as  opposed  to  a  four-class  study. 
These  SVM  limitations  are  here  explained  by  the  used  kernel 
drawbacks  as  the  considered  DTWK  is  not  positive  definite,  pre¬ 
senting  difficulties  in  its  applications. 

The  next  step  of  this  framework  is  to  explore  the  selection  of  the 
acceleration  criterion  to  have  the  best  compromise  between  model 
accuracy  and  number  of  patterns,  with  an  automatic  process.  The 
SVM  results  also  induce  a  possible  way  to  improve  algorithm  per¬ 
formances,  by  the  use  of  other  DTW  based  kernels.  Moreover,  the 
natural  extension  of  this  study  is  to  adapt  these  methodologies  to 
obtain  linear  estimations.  This  would  permit,  during  electric  uses, 
an  online  precise  continuous  estimation  of  the  battery  performance 
instead  of  a  level  class.  It  would  also  be  interesting  to  test  this 
methodology  with  a  larger  fleet  to  observe  the  different  battery 
behaviors  and  exploit  further  various  data  to  obtain  a  better  global 
accuracy. 

5.  Conclusion 

The  aim  of  this  study  is  to  use  an  innovative  methodology  built 
from  a  real  use  signals  dataset,  permitting  an  online  estimation  of 
the  battery  performance,  without  requirement  or  specification.  The 
dataset  is  formed  from  real  vehicle  uses  and  is  used  to  create  model. 
A  preprocessing  step  permits  to  extract  acceleration  patterns  in 
these  measured  signals,  under  specific  extractions  criteria.  Once 
these  battery  signals  patterns  extracted,  the  battery  behavior  is 
highlighted  to  be  altered  by  the  performances  levels.  These  pat¬ 
terns  are  then  used  as  a  training  dataset  in  the  battery  performance 
classification  approach. 

By  carefully  exploring  the  signals  alterations,  we  have  illustrated 
an  innovative  methodology  to  demonstrate  that  simple  pattern 
recognition  and  classification  techniques  permit  to  diagnose  bat¬ 
tery  performance.  This  pattern  analysis  methodology  is  also  able  to 
detect  battery  fault  detection  and  permits  the  identification  of  an 
abnormal  battery  behavior. 

The  presented  methodology  demonstrates  the  classification  of 
the  battery  performances  in  different  class  levels  with  an  inter¬ 
esting  confidence  level.  We  continue  to  find  new  methodologies 
and  new  pattern  recognition  algorithms  to  improve  presented  ac¬ 
curacy  results. 

However,  the  final  aim  is  to  obtain  a  continuous  estimation  of 
the  battery  performance  in  real-time.  Thus,  the  natural  objective  is 
to  extend  this  pattern  signals  investigation  methodology  to  accu¬ 
rately  quantify  a  battery  performance.  The  signals  pattern  analysis 
and  the  relationship  to  battery  performance  degradation  are  the 
main  focus  of  our  current  works.  Our  ultimate  goal  is  to  extend  this 
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methodology  and  these  results  to  accurately  estimate,  online,  a 
battery  performance  value,  only  from  the  signals  coming  from  a 
real-life  use. 
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